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Q_ | 1 Abstract 

Qh. The analytization procedure which allows one to remove nonphys- 

ical singularities of the QCD running coupling constant a s (q 2 ) in the 
infrared region is applied to standard as well as to iterative solutions 
of the two-loop renormalization group equation. Non-leading at large 
rS | momentum nonperturbative contributions in a s (q 2 ) are obtained in 

an explicit form. The coefficients of nonperturbative contributions 
expansions in inverse powers of the squared Euclidean momentum 
are calculated. For both cases considered there appear convergent at 
q 2 > A 2 power series of negative terms with different dependence on 
term numbers. 
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1 Introduction 

The nonperturbative contributions to observables as well as direct to the 
QCD running coupling constant a s have been widely discussed in recent 
years. The evidence provided by very different approaches suggests the exis- 
tence of power type nonperturbative corrections which reflect the complicated 
structure of the QCD vacuum. 

The nonperturbative contributions arise quite naturally when one uses an 
analytization procedure. The main purpose of this procedure is to remove 
nonphysical singularities from approximate (perturbative) expressions for the 
Green functions of QFT. The idea of the procedure goes back to Refs. |[], Q 
devoted to the ghost pole problem in QED. The foundation of the procedure 
is the principle of summation of information derived from the perturbation 
theory under the sign of the Kallen - Lehmann spectral integral. In recent 
papers [|], [§] it is suggested to solve the ghost pole problem in QCD by 
demanding the a s (q 2 ) to be analytical in q 2 (to compare with dispersive 
approach ||). As a result, instead of the one-loop expression a^'(q 2 ) = 
(4tt / bo) / \n(q 2 / A 2 ) which takes into account the leading logarithms and has 
the ghost pole at q 2 = A 2 (q 2 is the Euclidean momentum squared), one 
obtains the expression 



An 



1 A 2 



ln(g 2 /A 2 ) A 2 -q 1 



(1) 



Eq. ([!]) is an analytic function in the complex g 2 -plane with a cut along 
the negative real semiaxis. The pole of the perturbative running coupling 
at q 2 = A 2 is canceled by the nonperturbative contribution (A 2 | 9 2^ ~ 
/x 2 exp{ — (47r) 2 /(& fi' 2 )}) and the value 0^(0) = 47r/6 appeared finite and 
independent of A. The nonperturbative contribution in Eq. ([I]) can be pre- 
sented as convergent at q 2 > A 2 constant signs series in the inverse powers 
of the momentum squared. 

At the two-loop level one cannot succeed in finding the nonperturbative 
part of a^(g 2 ) in an explicit form. There were obtained the approximate 
formulas [f|, |6| which, in spite of their accuracy, do not provide an exhaustive 
information on the behavior of the nonperturbative contributions. The most 
important feature of the analytization procedure discovered is the stability 
property of the value of the analytized coupling constant at zero with respect 



to high corrections, «^(0) = «^(0) = 0^(0). This property provides a 
definite high corrections stability of a an (g 2 ) in the infrared region in a hall. 
In this paper we find out what happens with power corrections of non- 
perturbative ultraviolet "tail" of the analytized running coupling when going 
from the one-loop level to the two-loop level. 

2 The extraction of nonperturbative contri- 
butions 

Let us consider the two-loop /5-function 

M) = fog 4 + fog 6 , (2) 

where the coefficients 

1 2 

& = -74^4&i, &i = 102-yn/ (3) 

do not depend on the renormalization scheme choice. We can present the 
Gell-Mann - Low equation 

n 9 -^ = m, (4) 

u = q 2 /n 2 in a "rationalized" form by introducing the function a(x) of x = 
q 2 /A 2 of the form 

The two-loop differential equation for the running coupling constant can be 
written as follows: 

x^P- = -a 2 (x)-ba 3 (x), (5) 

where 

ft = 102-f„, 

Bl (11 -ln,f (b) 



For the case of three active quark flavors n/ = 3, 60 — 9, and b = 64/81 ~ 
0.7901. Integrating Eq. @ one obtains 

-61n (l + -—-]= lnx + lnc, (7) 



a(x) \ ba(x) 

where c being the integration constant. The transcendental Eq. (^) defines 
implicitly the function a(x). The real solution of Eq. (|7|) has two branches, 
the former corresponds to a positive decreasing at x — > 00 function a(x), 
whereas for the latter a(x) is negative and at x — > 00 it goes to —1/b. 
The positive decreasing at x — > 00 solution corresponds to the asymptotic 
freedom property, and this solution can be fixed further by the choice of the 
integration constant c, which means, in practice, the redefinition of A. The 
exact solution of Eq. (0) can be written |7j using the Lampert function W(y) 
defined by the equation y = W(y) exp{W(y)}. 

We shall consider the following two forms of approximate solutions of 
Eq. (|7|). The first one is the standard two- loop coupling 

a (2) (*) = y-^rnZ, (8) 



') 



where I = - In x, and the second one is the iterative solution || of the form 

a*(s) = - -X- /V (9) 



Z + Mn(l + 



This solution, when expanding in inverse powers of logarithms apart from 
the terms Eq. (^), gives rise to the term blnb/l 2 and additional terms of the 
form l~ 3 In 1 + ... which should be taken into account at the three- loop level. 
The exact solution of Eq. (|7j) also does not claim to be the description of the 
three-loop terms. At large x the functions (§), (|9D behave similarly but for 
small x the behavior is different and at x — 1 they have singularities of a 
different analytical structure. Namely, at x a 1 

a {2) (x)^-— b ]n(x-l), a^(x)^—±— (10) 

[x — \y 2(x — 1) 

These singularities are stronger than the integrable singularity of the exact 



solution of Eq. (0) of the form 1/yX — 1/c. But this is not an obstacle to 
the application of the analytization procedure. 



According to the definition the analytized running coupling is obtained 
from the initial running coupling by the integral representation 
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P(*) 



(11) 



where p(a) = lma an (—a — iO) = Ima(— a — iO). As initial expressions we will 
consider Eqs. (§), @. By making the analytic continuation of these equations 
into the Minkowski space x = —a — iO, one obtains, correspondingly 



o^(-cr-zO) 



a u (-a-iO) 



In a — in (In a — in) 2 
1 



In (In a — in) 



(12) 



In a - in + b In ( 1 + \ In a - f 



lncr + MnJ(l + i lncr) +^ 



7r + 6 arctan rrr — 



;i3) 



For the spectral function p(cr) it follows that 

7T 6 



P (2) (-) 



ln^ cr + 7r 2 
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(in 



2 cr + 7r 2 



X 



(in 2 



a — n ) arctan ■ 



71 



lncr 



2-7T In cr In v In a + n 2 



P» 



n + b arctan t-t — 



lncr + 61nJ(l + i lncr) + $ 



(14) 
(15) 



7T + b arctan 7-7— 



By the change of the variable of the form a = exp(t), the analytized expres- 
sions are derived from (|i~4|), ( |i~5"t ) as follows: 



a an( X ) = ~ dt 



■X 



00 



7T 



+ 



t 2 + Ti 2 (t 2 + 7l 2 Y 



(t 2 - 7T 2 ) arctan 2vrt In \Jt 2 + n 2 



(16) 



<n(z) 



dt 



it + b arctan -^ 



tt J x + e 

-oo 



t r 



t + b\n\ 1 + 



7T^ 

6- 



7T 



6 arctan 



6+i 



(17) 

Eqs. (0), ( |17D can be used to study a an (x) by numerical methods in par- 
ticular. We are interested in nonperturbative contributions which decrease 
at x — > oo much faster than the perturbative ones, therefore, their extrac- 
tion by a direct comparison with the initial expressions is difficult at large 
x. We shall obtain an explicit analytic formulae for the nonperturbative 
contributions as a difference between output and input expressions of the 
analytization procedure. 

Let us see what the singularities of the integrands of ([16]), ([H]) in the 
complex t-plane are. First of all, for both cases the integrands have the 
simple poles at t = lnx ± itt(1 + 2n), n = 0,1, 2, .... All the residues of the 
function exp(t)/(x + exp(t)) at these points are equal to unity. Making use 
of the formula 

(18) 



7T 

arctan — 
t 



li t 



m 
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one can obtain 



if 



7T' 



tc arctan ■ 



7T 



2vrf In Vt 2 + n 2 



ln(i — in) ln(t + in) 



(19) 



{t-iirf (t + in) 2 

This formula allows one to make transparent the singularities structure of 
the integrand in fll6"l) . We obtain 



a an \ X ) 
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2n~i 



dt- 



x + e* 



-b 



ITX 



m 



ln(£ — in) ln(t + in) 
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(20) 



For the second case under consideration, it follows from fll3|), (|17|) that 



O) 



1 

2n~i 



dt 



x + e* 



■x 



X 



(21) 



1 1 

t - itx + 6 In (l + \ - f ) t + in + b\n(l + l + f] 

The integrands in (J2TJ|), fl2"ip multiplied by t go to zero at | t |— > oo. This 
allows one to append the integration by the arch of the "infinite" radius 
without affecting the value of the integral. Closing the integration contour 
C in the upper half-plane of the complex variable t, by using the residue 



theorem for the one- loop part of ( fZCf ), one readily obtains the well-known 
result (|) 



offiC*) 



1 /• °° 

— - / dtf(t) = ^Res/(t = lnx + ivr(l + 2n)) + 



c n=0 



+ Res / (t = wr) = -L + -J— , (22) 

In x 1 — x 



where /(£) is the integrand of fl20|) at 6 = 0. The two- loop part of the 
integrand fl2"0| ) apart from the simple poles at t = lnx ± iir(l + 2n), n = 
0, 1,2, ... has logarithmic branch points at t — ±27r which coincide with the 
second order poles. Let us cut the complex t-plane in a standard way, t = 
±i7r — A, with A being the real parameter varying from to oo. For the 
case of coincidence of the branch point and the pole of the second order the 
following theorem can be proved. 

Theorem. If f(z) analytic function inside the circle of the radius r with 
the center at zero, then 

' dz ^f(z) = -2m I J ^ [f(x) - /(0) - z/'(0)] - 1/(0) - /'(0) In a 

C \ -a 

(23) 
where the contour C goes from z = —a — iO along the lower side of the cut 
till the circle of the radius 5 with the center at zero (r > a > 5 > 0), then 
goes by this circle around z = 0, and then along the upper side of the cut it 
goes to z = —a + iO. 

For the terms of (p0|) which are proportional to b, we also close the contour 
in the upper half-plane of the complex variable t excepting the singularity 
at t — Z7T. In this case an additional contribution emerges because of the 
integration along the sides of the cut and around the second order pole. The 
substitution z = t — m brings this contribution to the form corresponding to 
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the theorem fl2"3|) with f(z) = (b/2iri) exp(z)/(x — exp(z)). We shall call this 
contribution together with the one-loop contribution from the pole at t = iir 
as nonperturbative because the remaining contribution of the poles at t = 
lnx + 27r(l +2n), n — 0, 1, 2, ... results exactly in perturbative expression (§). 
In fact, according to the residue theorem 



«2*(*) 



a^ pt (x) 



-&£ 



?1=0 



ln(ln x + 2nin) ln(lnx + 2ni(n + 1) 



(lnx + 27rm) (lnx + 2m(n + 1)Y 



1 



lnx 



-b 



ln(lnx) 
(lnx) 2 



The nonperturbative part of the analytized coupling using Eqs. (|2 
can be derived with the result 



(24) 
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t 2 xe t 



x 



(25) 



Integrating by parts and making the change of variable a = exp(— t), we can 



represent (]25 ) in the form of one integral with finite limits 



a£> npt (x) 



Th +b \a^ +x I/ aH - lna) jx^w} 



(26) 



This formula is convenient for the numerical study whereas formula (|2 
suitable for finding the explicit power nonperturbative corrections at x 



) is 

oo. 

Let us now turn to the evaluation of the integral ([H|) for the analytized 
iterative coupling. The integrand has the simple poles at t = lnx±i7r(l + 2n), 
n = 0, 1,2, ... and at t = ±iir. It has also the logarithmic branch points at 
t — — b ± in. To provide the single-valuedness of the integrand, we draw two 
cuts beginning at t = — b ± iir and going to infinity to the left in parallel to 
the real axis of the complex t-plane. We close the integration contour by the 
"infinite" semicircle in the upper half-plane with going around the branch 
point by the sides of the cut and integrate along this contour. As in the 



previous case, all the terms in the sum of residues at t = lnx + m(l + 2n), 
n = 0, 1, 2, ... mutually cancel apart from the term of n = 0. Integrating also 
along the sides of the cut, we obtain 



«1(X> 



1 



+ 



lnx + MnM+ilnx) 2(1 - x) 



oo 

+ J ^ 1 - are 6 ^) (T 



e-hiO +vr 2 



(27) 



Performing the change of variable a = exp(— ^), we come to the integral with 
finite limits. For additional "nonperturbative" output of the analytization 
procedure we obtain the following representation 



Aa±(x 



1 f da 

+ 



2(1 - x 



a 1 - x{a/e)- b [1 - \n(-a\na)] 2 + vr 2 ' 



(28) 



3 Nonperturbative contributions at large q 2 . 

Consider the large x behavior of the nonperturbative contributions. It is 
seen from Eq. ( pip that this behavior is regular. Expanding the integrands 
in powers of 1/x and integrating, we find 



a^ np \x) 






n=l 



where 



-1 + bn 



n f dt lnte~ nt + l 



(29) 



(30) 



Performing the t integration in (|3t|) yields the following simple formula for the 
coefficients of the ultraviolet expansion of the nonperturbative contributions 

c n = — 1 + bn (1 — 7 — Inn) , (31) 

where 7 is the Euler constant, 7 ~ 0.5772. The leading term is 

ci = -1 + 6(1-7), (32) 



and the passing from the one-loop level to the two-loop one results in some 
compensation of the leading at large x term of the form 1/x. For the next 
terms there is no compensation, instead they increase with n increasing. 
Estimating at rif = 3 the first three terms of the series (p9[) , we obtain 



0.666 1.427 2.602 



a mn Pt{x) ^ _^ _—_ — __ (33) 



Let us consider the large x behavior of the "nonperturbative" contribution 
in (p7|). Expanding the integrand in powers of 1/x, we also obtain the infinite 
terms series 



oo it 



A <n(^) = E ^ (34) 

where 



n=\ x 



i — 1 

4* = -- - /d£e- bn(1+5) l . (35) 

Expression (|35| ) for the expansion coefficients is not so simple as (|3T| ) but it 
is explicit and can be calculated up to an arbitrary accuracy. We see from 
Eq. fl35|) that all c 1 ^ are also negative tending fast to —1/2 when n increases. 
Substituting the denominator in the integrand of (|3^) by its minimal value 
7r 2 + 4, one finds 

1 1 p~ bn 

At rif = 3 it gives for the first coefficient | df + 1/2 |< 0.0414. For the first 



three terms of the series (33), one can obtain 



A ,, . , 0.535 0.508 0.502 

Afl lW - "— -2 ^i- - - " (37) 

Thus, passing from the one-loop level of the analytization procedure to the 
two-loop level leads to approximately two times decrease of each term of the 
expansion (|3~4f) . 



4 Conclusion 

For both cases considered the nonperturbative contributions can be presented 
in the form of series in inverse powers of the momentum squared, all the terms 

9 



are negative. Although the convergence radii of the series (p9[), ( |34"D are the 
same and equal to unity, the dependence of the expansion coefficients on 
term numbers is different. This is defined by the expressions taken as an 
input for the analytization procedure. Choosing as an input the standard 
two- loop expression (f|), we obtain the "standard nonperturbative correc- 
tions" of the form (^). For sufficiently large x = q 2 /A 2 the nonperturbative 
contributions are defined mainly by the leading terms in (|33|) , (|37|). The 
partial compensation of the leading terms found for both cases considered 
when we pass from the one-loop level to the two-loop level can point to the 
tendency of a high-loop minimization of the nonperturbative contributions 
in the ultraviolet region. 

The analytization procedure may be not an ultimate step of construction 
the "physical" running coupling constant. In terms of Ref. || the analytized 
coupling constant corresponds to the "regularized perturbative part" of the 
full coupling constant which contains also the "genuine nonperturbative" 
contribution. As a step in the direction of the full running coupling con- 
stant it has been suggested in |], [UJ that the one-loop analytization output 
be modified in a minimal way by introducing the two additional pole type 
nonperturbative terms to provide the ultraviolet convergence of the gluon 
condensate. In this case the model for the running coupling constant arises 
with the enhancement at zero momentum and dynamical gluon mass m g . 
This mass can be fixed |~l| by the condition of minimum of the nonpertur- 



bative vacuum energy and is estimated m g ~ 0.6 GeV for the "standard" 
value of the gluon condensate (0.33 GeV) 4 . The next step of the model 
building for the full running coupling constant on the base of the two-loop 
analytized running coupling is under study. 

I am indebted to B.A. Arbuzov, V.A. Petrov, V.E. Rochev for useful 
discussions. This work was supported in part by RFBR under Grant No. 99- 
01-00091. 
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